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I. INTRODUCTION 

Interferometry plays a central role in physics, with ap- 
plications ranging from sensitive phase measurements, 
e.g. to monitor the spatial displacements experienced by 
gravitational wave detectors [T], to fundamental tests of 
quantum physics, such as the wave-particle duality of 
increasingly complex quantum objects [21 [3] . Given the 
power of interference experiments, it is natural to ask how 
their scope can be extended to access entanglement 
the second pillar of nonclassicality [1], and an impor- 
tant resource in quantum information [HHI] • Such a com- 
bined witnessing of spatial interference and entanglement 
would not only amount to a striking demonstration of 
the departure of non-local quantum behavior from clas- 
sical physics. It might also be used for entanglement- 
enhanced metrological applications such as phase esti- 
mation schemes or quantum lithography [5J [9] ■ 

So far all experiments that combined entanglement 
with spatial interference were based on photons [5lU2|. 
This is due to their great practical use in information 
science, and above all due to the existence of a mature 
technology for producing and manipulating photonic sys- 
tems. However, the recent advances in the control of ul- 
tracold atoms suggest that it will be possible to carry out 
similar experiments with material particles. In particu- 
lar, tailored Einstein-Podolsky-Rosen (EPR) entangled 
atom pairs can be produced by dissociating Feshbach 
molecules [P3HI6] or by colliding Bose-Einstein conden- 
sates P~7l420| . in a process similar to parametric down 
conversion of laser light, the established method to gen- 
erate entanglement among photons [2T] . 

It would be a great experimental advancement to 
demonstrate a nonlocal spatial interference effect with 
particles of matter. This would establish the presence 
of both entanglement and the wave-particle duality in 
a single experiment on tangible material objects, allow- 
ing one to transfer potential application schemes such as 
quantum lithography [5J [S] from photons to the realm 
of quantum matter. Moreover, interpretational issues 
such as the nonlocality of Bohmian trajectories could be 
addressed [221 • However, it is not obvious whether the 



observation of an interference pattern in the coincidence 
signal of two detectors already proves the existence of en- 
tanglement. This is the more so, as any experiment will 
be characterized by a non-ideal EPR source and other 
imperfections leading to a reduced fringe visibility of the 
interference signal. 

In this article we provide an entanglement criterion 
for nonlocal spatial interference based on EPR entangled 
particle pairs, and we work out the conditions for the 
successful detection of a nonlocal interference signal. A 
schematic generalizing the single-particle Young experi- 
ment is depicted in Fig. [T] A source creates a pair of 
EPR entangled particles which travel freely into opposite 
directions until each one passes through a grating or dou- 
ble slit. After a further free evolution their positions are 
recorded by spatially resolving detectors located at oppo- 
site sides of the experiment. Even for a source emitting 
ideal EPR particle pairs no interference will be observed 
at each of the detectors. However, a 'nonlocal' interfer- 
ence pattern is expected to emerge if one analyzes the 
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FIG. 1. Generic setup of a nonlocal Young experiment with 
EPR entangled particle pairs. A source provides EPR en- 
tangled particle pairs of equal mass, which are then sepa- 
rately passed through a grating or double-slit (with slit sep- 
aration d and slit width a). While the position measure- 
ments yield an unstructured pattern on each screen, inter- 
ference fringes are expected if one evaluates the correlated 
outcomes x± + X2- Taking the EPR correlation to be in the 
transversal x-direction it is sufficient to describe the longi- 
tudinal 2-motion, which separates the particles, by classical 
counter- propagating trajectories. 
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combined detection records at both sides by focusing on 
the center-of-mass of the coincident pairs. 

A major challenge in such nonlocal interference experi- 
ments lies in verifying the presence of continuous variable 
entanglement. Once the gratings are passed the two- 
particle state of motion is strongly non-Gaussian, such 
that the standard Gaussian entanglement criteria [2"5H2"5] 
cannot be used. Moreover, only position measurements 
are easily doable with material particles, practically rul- 
ing out tomographic techniques of entanglement verifi- 
cation. Recently, we described a viable method based 
on modular variables, which captures situations where 
entanglement gives rise to spatial interference |26| . In 
the following we introduce the corresponding criterion, 
adapted to the specific correlations displayed by the EPR 
interference experiments. 

The article is structured as follows: In Sect. [II] we de- 
scribe how the nonlocal interference pattern can be cal- 
culated for finite EPR sources, allowing us to discuss the 
requirements and conditions for observing nonlocal inter- 
ference. After briefly explaining the concept of modular 
variables, the entanglement criterion is then formulated 
in Se ct. |III| along with a discussion of its robustness. In 
Sect. |IV| we discuss different experimental scenarios, be- 
fore presenting our conclusions in Sect. [Vj 



II. NONLOCAL SPATIAL INTERFERENCE 
FROM EINSTEIN-PODOLSKY- ROSEN 
CORRELATIONS 
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FIG. 2. Correlations expressed by the squeezed Gaussian |T]) 
representing a normalizable EPR state. 
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The resulting correlations are sketched in Figure [2] for 
both <7 X! rel an d CTp,cm much smaller than the uncertainties 
of an (unsqueezed) minimum uncertainty state. 



A. The normalized EPR state 

Einstein, Podolsky, and Rosen considered originally 
the idealized state [I] 



dar|a;)i|x)2 = / dp\p)i\ - p) 2 



Switching to center-of-mass coordinates, and disregard- 
ing normalization, it may be written as 
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The state supports perfect correlations both in the rela- 
tive position x rc i = xi — X2 and in the center-of-mass mo- 
mentum p cm = p! + p 2 . The expectation values of these 
commuting observables vanish for the idealized EPR 
state, as do their variances a 2 x rel and Cp Cm . This implies 
that the conjugate relative momentum p re i = (pi — p 2 )/2 
and the center-of-mass position x cm = (xi +x 2 )/2 remain 
undetermined. 

The idealized EPR state is readily generalized to a nor- 
malized squeezed Gaussian wavefunction exhibiting finite 



B. Interferometric two-particle evolution 

We proceed to determine the joint probability distri- 
bution for recording the two EPR particles on detection 
screens after each passed a grating with N slits. To this 
end it is helpful to divide the evolution into three steps: 
(i) the free time evolution from the source to the grat- 
ings, (ii) the instantaneous effect of passing through the 
gratings, and (iii) a further period of free time evolu- 
tion from the gratings to the screens. This can be done 
since all the correlations probed in such a setting reside 
in the transversal motion of the two particles. In the 
Fraunhofer approximation the longitudinal motion from 
the source to the screens may be viewed as taking place 
with definite velocities which are to be averaged in the 
end. The longitudinal position can thus be considered as 
a parametrization of time. 

We take the longitudinal direction to be the z-axis and 
the grating bars to be aligned along the y-axis such that 
the relevant motion takes place in the x-direction. The 
free time evolution of the Gaussian wavefunction (JfJ from 
the source to the gratings can be determined analytically. 
Denoting the arrival time of the particles at the gratings 
as T, the evolved state immediately before passing the 
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gratings reads as 




d#cmdx re l |#cm)cm|^rel)rel 
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Here we transformed from the center-of-mass momentum 
Pcm to the center-of-mass position x cm , with a x . cm = 
hj (2(T p . cm ) the corresponding (large) uncertainty. The 
evolution introduces additional phase factors and results 
in a dispersion-induced broadening of the original Gaus- 
sian wave packets described by the complex dispersion 
factors 



1 



iht 



4 < cm m' 



6,rcl — 1 + 



iht 



'x,rcl' 



(3) 
(4) 



associated with total mass 2m and reduced mass m/2, re- 
spectively. The normalization factor Mr and the global 
phase ipT, which are readily calculated, will not be re- 
quired in the following. 



The effect of traversing the gratings can be captured 
by the grating operators 



dx\nd + x)i(nd + x\ 



(5) 



which are projectors acting on the individual particles, 
i = 1,2. The gratings are taken to consist of N slits, 
with slit distance d and slit width a, see Fig. [T] The slit 
indices are taken from 
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which guarantees that the gratings are arranged symmet- 
rically with respect to the z-axis for all N. The gratings 
are assumed to be ideal in the sense that imperfections 
and the dispersion force between particle and grating can 
be neglected. 

Immediately after traversing the gratings the state fol- 
lows from the projection Gi ® G 2 |'I , t)- Switching to 
the particle coordinates x\ = x cin + x rc \/2 and X2 — 
x C m — a^rci/2 the entangled two-particle thus takes the 
form 
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C. Conditions for nonlocal interference 

We can now identify the conditions for observing non- 
local spatial interference at the screens. As follows from 
the general discussion of two-particle correlations in [55] , 
nonlocal spatial interference requires that the state pre- 
pared by the gratings is correlated with respect to the 
slits traversed. That is, if we measure which slit a parti- 
cle took on one side, by detecting it right after the grat- 
ing, we must be able to infer which slit the other particle 
took on the other side. This is guaranteed by requiring 
that 

OVrcl|£T,rol| < d, (8) 

because the second Gaussian in Eq. ^ can then be well 
approximated by a Kronecker delta for the slit indices 
n and n' (along with with a factor describing irrelevant 
intra-slit correlations) . 

The slit correlation condition Q combines two require- 
ments. On the one hand, the initial state must be suffi- 
ciently squeezed in the relative position, a x ro i <C d. On 
the other hand, to guarantee (fsj) at the relevant time T 



when the gratings are passed, the dispersive broadening 
of the wave packet evolving from the source to the screen 
must remain sufficiently bounded. The propagation time 
T from the source to the gratings thus cannot exceed 

T max = tT 2 x rcl m/h < d 2 m/h. (9) 

This shows that there is a tradeoff between the squeezing 
of the relative position a x ,rei and the maximum admissi- 
ble propagation time from the source to the screen. This 
is no problem of principle since the longitudinal motion 
might be uncoupled from the transverse dynamics, allow- 
ing one to propagate the particles with arbitrary velocity 
to their gratings. In practice, one may well be forced to 
start out with an isotropic two-particle state, using aper- 
tures for defining the longitudinal direction. In this case 
the dispersion tradeoff ^ can be a quite severe restric- 
tion. 

As a second condition one must require that the N slits 
are illuminated uniformly by the incoming wavefunction. 
This is ensured by the first Gaussian in Eq. |7]) provided 

Ifr.cml > Nd. (10) 
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Since the slit correlation condition ([8]) already requires 
that the relative dispersion remains modest, T < T max , 
we can estimate 



because a x rc i <C cr x 




(11) 



The initial state (fll) must there- 



fore exhibit a large center-of-mass uncertainty cr x cm > 
Nd. 



D. The modular momentum entangled state 



Applying the conditions Q and ( 10 ) to Eq. Q yields 
a greatly simplified expression for the slit-entangled two- 
particle state present once both gratings have been 
passed, 
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It describes a superposition of the particle pair pass- 
ing through opposite slits. This becomes clear once we 
rewrite the state in position representation, 



Xi,x 2 \^' T ) — -j= (xi - nd, x 2 - nd\^ s 



(13) 



where the interference pattern is captured by the fringe 
function 



N-l 



F N (^) = 1 + -^2(N- j) cos(27rjO- 



(17) 



It reduces to a sinusoidal fringe pattern in the case of 
double-slits, while F^(x) develops sharpened main max- 
ima and suppressed side maxima for N > 2. Note that 
the period of the fringe pattern is given by the 'grating 
momentum' h/d. 

A distinct interference pattern can only emerge if the 
envelope in (16), as given by the momentum distribu- 
tion | {pi , jc?2 1 ^T) | 2 j varies slowly over the extension of a 
single period h/d, and if it is sufficiently broad to cover 
several fringes. These conditions are met in the present 
case since the width of the momentum distribution of 



( 14 ) is essentially determined by single-slit diffraction, 



i. e. by h/a. This is always greater than the grating mo- 
mentum h/d, since d > a. More precisely, the envelope 
is determined by a convolution of the opposite slit pair 
contributions and the correlated relative motion, 



(pi,p 2 \* s ) = 




(18) 



xcxp 



Here the normalized wavefunction |\E' S ) describes that 
both particles are confined to a single pair of opposite 
slits, 



(xi,x 2 \ i &s 



exp 
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with Xa( x ) = © (f + x ) © (f — x )- As before, the phase 
ip s and the normalization factor Af s will not be required 
in the following. 



The wavefunction ( 13 ) is an instance of a modular mo- 



mentum entangled (MME) state, a class of states dis- 
cussed in a more general context in |26| . As for any 



MME state, the momentum representation of ( 13 1 reads 



ipt,P»\%.) = -L Z-' i{Pl+P * )nd/ \Pl,P2\^s). (15) 

This in turn implies a nonlocal interference behavior if 
the particle momenta arc measured. The joint momen- 
tum probability distribution takes the form 

\(puP 2 \%)\ 2 = \<Pi,P2\*.)\*F N ( {Pl + J 2)d ) , (16) 



E. Far-field interference pattern 

The discussed momentum interference effect is easily 
observed by letting the particles propagate freely for a 
time T 2 from the gratings to remote detection screens 
on each side. Denoting the final state of the particles 
as |W /), the joint spatial detection probability is directly 
determined by the momentum distribution (16) of \^' T ), 
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This assumes that the screens are placed sufficiently far 
away from the gratings such that one is in the dispersion- 
dominated limit, T 2 3> mN 2 d 2 /h. 

The position measurements at the detection screens 
may thus be viewed as effective momentum measure- 
ments on \^' T ). It follows that the joint spatial prob- 
ability distribution reproduces the nonlocal momentum 



interference pattern (16), 
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(20) 
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This result exhibits the expected nonlocal interference 
behavior. No fringe pattern will be visible if one looks 
at either of the screens since the integration over the 
unobserved particle position will remove the fringe func- 
tion in (161, leaving only the broad envelope determined 
by |^ s ). Only by recording the coincident detections at 
both screens and by collecting the center-of-mass posi- 
tions x\ + x 2 will an interference pattern emerge. 

This proves that it is possible to establish nonlocal 
interference by exposing EPR entangled particle pairs to 
gratings, and in this sense to perform an entangled Young 
experiment. The expected spacing of the fringe pattern 
is given by T 2 h/md. 

We identified the slit correlation condition (fsl) and the 



with 



uniform slit illumination ( 10 1 as requirements for a suc- 



cessful implementation. In the following section, we will 
show how the nonlocal spatial interference pattern (20) 



can serve as the basis for a rigorous verification of the 
underlying entanglement. 



III. INTERFEROMETRIC ENTANGLEMENT 
VERIFICATION 

As impressive as the correlations expressed by the non- 
local interference pattern (20) may be, it is not clear a 



priori that they cannot just as well emerge from a classi- 
cally correlated quantum state, without resorting to en- 
tanglement. To exclude this possibility, one must testify 
the presence of entanglement with a suitable entangle- 
ment criterion. Ideally, this should not require measure- 
ments beyond the ordinary position measurements giving 
rise to the interferometric correlations of the EPR Young 
experiment; in particular we should avoid the necessity 
of an unfeasible continuous variable state tomography. 



A. Modular variables. 

Such an entanglement verification can be achieved us- 
ing an entanglement criterion in terms of modular vari- 
ables The latter proof useful to capture spatial in- 
terference phenomena [2TH3T3] . They formally decompose 
the position and momentum operators into step-like in- 
teger components and N p and sawtooth-like modular 
components x and p, 



x=N x d + x, p=N p - + p. 
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(21) 



Expressed in the position eigenbasis the modular and the 
integer position are thus given by 



da; x (x) \ x) (x |, 
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The momentum operators are defined similarly by the 
spectral function 
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Note that the standard position and momentum eigen- 
vectors |a;) and \p) can also be interpreted as the joint 
eigenstates of the respective integer and modular observ- 
able. The latter can thus be deduced from ordinary po- 
sition and momentum measurements. 



B. Moments and variances 

Interpreted in terms of the modular variables, the cor- 
relations displayed by the MME state ( 13 1 describe re- 



duced variances for the total modular momentum p cm = 
Pi + P~2- This is the reason for naming it modular mo- 
mentum entangled. Also the spread of the relative integer 
position N x rc i = N Xi i — N x 2 is reduced. In particular, 
if the slit correlation condition ^ is satisfied the vari- 
ance of N^ rci vanishes by construction, as a result of the 
opposite slit pair correlations, 



((*T|N,, r cl|*T» 5 



= 0. 



(26) 



Using (13), (14) and noting N x (x) = (x — x{x))/d one 



finds that in fact all moments m > 1 vanish, 
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On the other hand, in the relevant limit a <C d and using 
(16), the moments of the total modular momentum are 



given by 
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Here we used the scale separation between the width of 
the fringe pattern envelope and its period; this permits 
one to apply the approximation 
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which is valid for a A-periodic function x( x ) an( i an en ~ 
velope function E(x) that varies slowly over the extent 
of a single period A. 



Putting the fringe function (17 1 into (28 1, one obtains 
for the first two moments 
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The positive function 5*2, also found in 
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is defined by 
(32) 



3=1 



It is bounded by 1 > ^(AT), increases monotonically, 
and is well approximated by its asymptotic form 
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involving Euler's constant 7 c 
the total modular momentum 



0.577. The variance of 
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(1 - S 2 (N)), 



(34) 



thus decreases with a growing number N of slits. For 
large N this squeezing of the total modular momentum 
scales as 



<(Ap cm ) 2 
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C. Post-measurement analysis 

In the proposed Young test the experimenter perfoms 
ordinary position measurements directly behind the grat- 
ings and in the far field, yielding the joint probability 
densities prob(xi, x 2 ) and prob(pi,p 2 ), respectively. The 
required modular variances ((AN^.td) 2 ) and ((Ap cm ) 2 ) 
are then obtained by a post-measurement analysis of this 
data. 

Specifically, the moments of the relative integer posi- 
tion are obtained by evaluating 



(N™ rel )=/ d Xl dx 2 [N x ( Xl ) ~ N x (x 2 )] m 
J— 00 J —00 

xprob(a;i, x 2 ) (36) 
with N x (x) = (x — x{x))/d. Similarly, one calculates 
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where p(p) is defined in (25). From an operational point 



of view, this is all that is required to test the entangle- 
ment criterion (41). 



D. Shifted modular variables 

Carrying out the post-measurement analysis one can 
always choose the position and momentum coordinates 
in such a way that the maxima of the interference pat- 
tern coincide with vanishing values of the corresponding 
modular variable. This reflects the optimal choice and 
was the case in our calculations so far. 

The general case can be modeled by introducing an 
additional phase ip into ( 20 ) which shifts the interference 
pattern. The expression ( 28 ) for the moments of the total 



modular momentum then becomes 
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This results in the modified variance 
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(39) 



(40) 



A finite ip can result in a substantial deterioration of the 
total modular momentum squeezing, while the moments 
of the relative integer positions remain unaffected. In the 
remainder we consider again the case ip = 0. This is no 
restriction of generality due to the freedom of choice of 
ip in the post-measurement analysis. 



E. Modular entanglement criterion 

The reduced fluctuations in the integer relative posi- 
tion and the total modular momentum can be used to 
verify unambiguously the underlying entanglement. This 
is achieved with an entanglement criterion similar to the 
modular entanglement criterion derived in |26j . where 
the squeezing was considered to occur in a different set 
of two-particle observables, namely the modular relative 
position x Te \ = x\—x 2 and the total integer momentum 
N PiCm = N Ptl +N p>2 . 

In the present case the relevant entanglement criterion 
reads 



rf 2 
h 2 



((Ap ^ 



J cmJ IP 



((AN :c , rcl ) 2 ) p <2Cp, Nx . (41) 



Any state satisfying this condition must be entangled. 
Note that the criterion (41 ) is sufficient but not necessary. 
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The constant Cp,i\i^ is given by the smallest root /j,q of the 
equation 



_d_ 

d.r 
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(42) 



2=1/2 



with M(a, b\x) the Rummer function. Numerical evalu- 
ation yields Cp, Nj; 01)78 235. 



The MME state ( 13 1 satisfies the entanglement crite- 



rion (41) for any N > 2, as follows directly from the 



variances (26) and (34) 



(1 - S 2 (N)) + < 2C p m x 



(43) 



This proves that it is indeed possible to detect unam- 
biguously entanglement based on the nonlocal interfer- 
ence which is produced by exposing EPR entangled par- 
ticle pairs to a Young-like grating setup. Already for the 
least sensitive entanglement scheme, the case N = 2 of a 



double slit on each side, the squeezing function ( 32 ) eval 



uates as 5*2(2) = 0.30, resulting in a sum of uncertainties 
staying 25% below the threshold. 



F. Robustness of the entanglement detection 

To get an understanding for the robustness of the en- 
tanglement detection scheme one may ask how many clas- 
sical correlations can be admixed to the EPR state with- 



out compromising the criterion (41). To this end we in- 



troduce the classically slit-correlated state 
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Compared to the MME state pmme = \^'t)(^'t\ deter- 
mined by (13), the state (44) lacks the coherences be- 



tween different opposite-slit pairs. It does therefore not 
exhibit nonlocal interference. It carries the variances 
((AN Xircl ) 2 ) cl = and ((Ap cm ) 2 ) c i = h 2 /6d 2 ; this latter 
variance of the total modular momentum is the maximum 
possible, reflecting complete ignorance. 
Let us now consider the mixture 



(1 - MljpMME + Wp c i 



(45) 



of the MME state ( 13 1 and the classically slit-correlated 



state MM with to 6(0,1). We find that the variances 



evaluate as ((AN x rc 



and 



((Ap cm ) 2 )„, = (1 - W )((Ap cm ) 2 ) MM E + w((Ap cm ) 2 ) cl 
h 2 



6d 2 



[l-(l-w)S 2 (N)}, 



(46) 



noting that all involved first moments vanish, 

(N Xt rel)MME = (Nx.rel)cl HPcieImME = (Pcm)cl = 0. 



Comparing (34) and (46) one sees that the squeezing 



of the total modular momentum is diminished by the 



amount of classical admixture w. This corresponds to 



a reduced visibility, and the fringe pattern ( 16 1 gets re- 
placed by 



|<Pi,J>a|¥* 



w + (l-w) F J 



N 



(Pi + P2)d 

h 



(47) 



In the case of double slits, N = 2, we thus find that 



the entanglement criterion ( 41 ) remains satisfied as long 
as to < (12Cp,N K - l)/S 2 (2) + 1 = 0.79. 

In other words, we can admix up to 79% of a classi- 
cally correlated state and still detect the entanglement 
in the blurred fringe pattern. This robustness manifests 
the power of the modular entanglement detection scheme 
and it provides a comfortable cushion to deal with po- 
tential noise sources and experimental limitations, such 
as decoherence and a finite detection resolution, which 
reduce the fringe visibility. 

In the opposite case, where the source produces uncor- 
rected particle pairs, the wavefunction behind the grat- 
ings is described by the separable state 



|*scp) = |-0ms)l|-0m S )2 

with single-particle multi-slit states 

{x\ip ms ) = — L= ^ (x - nd\ip s 



(48) 



(49) 



Here the state \ip s ) corresponds to the single-particle 
state prepared by a single slit. The state l^scp) then 
leads to local interference patterns on each side. This im- 
plies correlations in the modular total momentum which 
reduce its variance. However, the lacking slit correlations 
result in a substantial variance of the relative integer po- 
sition 



(A 2 N a ,, rol ) sop = -(N 2 



1): 



(50) 



As 



in contrast to the vanishing variance ( 26 1 . 
pects, already for N = 2 this exceeds substantially the 
threshold value 2CpM x of the entanglement criterion (41 ). 



G. Imperfect EPR sources 

Another possible reason for a reduced interference vis- 
ibility are imperfect EPR sources. We therefore dis- 
cuss in the following how the conditions and results de- 
rived above are affected if the initial state is not a pure 
EPR state ([T]) but a mixture of such states. They will 
be characterized by the phase space coordinates T = 

(aJcm, #rei iPcm,Prci) indicating where each EPR state is 
initially located with respect to the center-of-mass and 
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relative coordinates. 



I* 



dx cm d.37 re l |#cm) cm |*^rcl) rcl 



(r) v _ 1 

EPR/ _ fey 

w Z7T G x ,0111^ x ,rel 

(^cm ^cm) (*^rcl *^rel ) 



xcxp 



An 2 An 2 



( (0) (0) 

x exp 1 V i-^r 



(51) 



with dr = dicmda;|°|dpcmdpj:"| . It is thus determined by 
the probability distribution function /x(r), taken in the 
following to be a Gaussian centered at the origin, which 

is fully characterized by the standard deviations cri°cm, 
rr {0) rr {0) and rr (0) 



Comparison with Eq. ([!]) shows that the previously con- 
sidered EPR state is centered at T = 0. The general 
mixture is given by 



p, = |dr M (r)|4 r 4)(4 r p ) R |. 



(52) 



To see the effect of the mixing ( 52 ) , we first determine 



the interference pattern of a (moderately) displaced, pure 
EPR state (51). Freely propagating the wavefunction 

I^epr) f° r time T and then through the gratings yields 



I* 



'(r)\ 



— > dx dx' \nd + x)i\n'd + x')2 exp i < 

JM' T ^-1 J-i J-% \ 



ft ~~\~ n> . x ~\~ x \ . f - , 

a H — | + i (p IC \ ( [n — n \ a + x — x ) 



xexp 



([n + n']d + as + x' - 2x { c °2 - P^T/m) 2 ([n - n']d + x-x'- xfj - 2p ( ^T/mf 



16 <cml6r,cm| 2 



4a 2 , rol |^ T ,rel| 2 



(53) 



Here, we introduced the phase functions for the center- 
of-mass and relative motion, 



<Acm(a?) = 



0rol(x) 



(0) 
Pern X 

^l£T,cm| 

(0) 
Prcl g 



%T(x 



r (0 h 2 



tSmo-icm 



hT(x 



l£r,cm| 
.(0)^2 



"rel - 



4m< rel |fr ire i| 2 - 



(54) 
(55) 



Comparing the form ( 53 ) of the displaced wavefunction 
with Q it follows that the slit correlation condition ^ 
and the requirement of uniform illumination ( 10 ) remain 
necessary conditions. 

In the following, we determine what additional con- 
straints must be satisfied by the displacements T for a 
successful Young test. First, it must be guaranteed that 
the grating is still illuminated uniformly. Noting that 
the wavefunction (53) is centered around the classical 
displacements 



r (T) 



,(0) 



T cn _ _(o) 

A rcl ~ X rcl 



Jo) T 

Pcm-L 

2m 

2 Prel± 

m 



we get the requirement 



<iVcL 



(56) 
(57) 

(58) 



Next, consider the impact of the displacements on the 
resulting slit correlations. To this end, it is helpful to 

(T) 



express sr^V in modular variables, x^ T ' 



"rcl 



Similarly to the undisplaced case, the first Gaussian in 
(53), combined with the slit correlation condition (Jsj) , 



implies ideal correlations, n — n' = Nf^L. Note that in 
contrast to the case T — 0, we must now also take into 
account that it is not necessarily opposite slit pairs that 
are correlated. Uniform illumination of the slits on both 
sides demands that the offset N^Xd is small compared 
to the extension of the grating Nd, 



nZI « N - 



(59) 



Moreover, to make sure that both particles can pass the 
gratings in spite of their correlation the modular part 
must satisfy 



AT) 



a 

< 2' 



(60) 



If the conditions ( 58 )— ( 60 ) are met, we obtain the in- 
terference pattern 



|( Pl ,p 2 |< r) )| 2 H(pi,p 2 |*s}| 

x F N 



(61) 



fal +P2\d _ Pern d 



(T) 

with an interference order N = N — \N^ rel \. It implies 
that the interference is remarkably robust against phase 
space displacements, since only a shift in the center-of- 
mass momentum, pim, directly affects the phase of the 
nonlocal interference pattern. 

With this we are in a position to discuss the implica- 
tions for the mixed state (52). The center-of-mass re- 



■'rel 



quirement (|58| leads to the constraints cri° C m *C Nd and 



Cp.lm <C mNd/T, and the condition (591 for the rela- a) 



tive motion implies <J^ le 

While these "classical" requirements are relatively easy 
to meet, the sensitivity of the interference pattern ( plj ) 
with respect to phase averaging demands a significantly 

tightened control of 

Specifically, the blurred interference pattern due to the 
phase averaging results in an increased variance of the 
total modular momentum as compared to (31 1, 



! < Nd and ojjj < mNd/T. 



(Pcm)^ 



6d 2 



1 



exp 



T,cm I 



S 2 (N) 



.(62) 



A significant reduction of the fringe visibility is thus to 
be expected once the the total momentum spread a p °l m 
exceeds the grating momentum h/d. This indicates the 
level of control of the initial state required for a successful 
entanglement detection. 



IV. EXPERIMENTAL IMPLEMENTATIONS 

Let us now turn to possible experimental demonstra- 
tions of nonlocal Young tests and of the verification of en- 
tanglement based on EPR correlated particle pairs. We 
will see that conceivable implementations are quite di- 
verse, which is a result of the generic nature of the pre- 
sented scheme. 



A. Photon experiments 




b) 




B{t) 

A>A> 



FIG. 3. Schematic setups for realizations of an entangled 
Young experiment based on the controlled dissociation of ul- 
tracold Feshbach molecules. An ultracold Feshbach molecule 
[35] initially trapped by two crossing laser beams is disso- 
ciated by a magnetic field pulse, which produces an EPR- 
entangled particle pair. In (a) the laser guide is switched 
off after the completion of the dissociation so that the atoms 
fall freely towards the gratings on each side. In (b) a se- 
quence of two magnetic field dissociation pulses generates a 
dissociation-time entangled particle pair [151 116] . Once the 
early and the late wave packets overlap due to the dispersive 
time-evolution, a nonlocal interference pattern can emerge in 
the recorded positions of the two particles. 



double slits, based on the criterion (41 1 obtained from 
ESI. 



While this article focuses on entangled matter waves, 
it should be emphasized that nonlocal interference can 
be observed with photons as well. Here one can rely on 
established methods for generating EPR entangled pho- 
ton pairs, e.g. by parametric down conversion of a laser 
beam |31l I32| . Such bipartite interference experiments 
have been performed in different contexts, for instance 
in quantum lithography 9 , ghost interferometry [33], or 
the spatial implementation of qubits 10J. These exper- 
iments demonstrated nonlocal spatial interference, but 
they could not verify the continuous variable entangle- 
ment for want of a rigorous criterion. 

The analysis in Sects. [TT] and |III| can be carried over 
to EPR entangled photon pairs because the interfero- 
metric arrangement allows one to treat the photons as 
distinguishable particles, to be recorded at distinct posi- 
tions Xi, x 2 on a photodetector. Moreover, the Kirchhoff 
diffraction integral in Fraunhofer approximation for the 
photonic modes yields the same expressions as the free 
propagator of a quantum particle in paraxial approxima- 
tion. All that needs to be done is to express the evolved 
time t in terms of the longitudinal momentum p z = zm/t, 
which is given for the photons by p z = h/X. 

In a recent article [34] Carvalho et al. describe an ex- 
periment with down-converted entangled photons. They 
report a significant observation of entanglement with 



B. EPR pairs from atomic BECs 

A recent proposal by Kofler et al. [20] sets out to pro- 
duce EPR entangled atom pairs of metastable helium 
atoms. This is based on a four-wave mixing processes in 
a Bose- Einstein condensate as in [36] . The helium atoms 
are kicked against each other by stimulated Raman tran- 
sitions and collide by s-wave scattering. If the applied 
laser pulses are sufficiently weak, such that on average 
only a single pair of atoms is detected in the end, the 
resulting two-particle state is well described by Eq. ([!]). 
The pair then falls freely under gravity, each of the atoms 
traversing a double slit aperture, until they hit the de- 
tector where they are recorded with high efficiency and 
resolution. 

As in most of the entanglement tests based on the cri- 



terion (41 ), the detector should be movable, to be placed 



alternately either directly behind the slits or sufficiently 
far away. When positioned close to the slits it serves to 
ascertain that the particle pair is sufficiently correlated 
with respect to the slits traversed. This is quantified by 
calculating the variance of N X)re i = N^i — from the 
observed positions. When positioned in the far field, the 
detector records essentially the transverse momenta p 1; 
P2 of the particles. After correlating the data of both 
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particles the resulting nonlocal interference pattern ( 19 ) 



will have a finite contrast which limits how well-defined 
the phase of the pattern is. This phase uncertainty is 
quantified by the variance of the modular part of the to- 
tal momentum p cm = p 1 + p 2 , the second ingredient to 
the entanglement criterion. 



C. EPR pairs from molecular Feshbach dissociation 

Another possibility to generate clouds of EPR entan- 
gled atom pairs is to use the controlled dissociation of 
Bose- Einstein condensed Feshbach molecules [LTHTol 155] . 
Starting from a sufficiently dilute condensate and apply- 
ing weak dissociation pulses, it is again possible to fo- 
cus on single EPR atom pairs by removing multiple pair 
events in a post-selection procedure. A detailed investi- 
gation of this dissociation scheme, including the confin- 
ing geometry induced by trap and guiding lasers, can be 
found in [THl [37] . Using the well-developed techniques of 
atom interferometry [2] one would then have to imple- 
ment the required gratings by material or light-induced 
structures. A schematic of a possible setup is given in 
Fig. [3] (a) . The horizontal propagation is induced by the 
dissociation process, while gravitation causes a vertical 
acceleration towards the gratings. The relevant trans- 
verse EPR correlations thus reside in the horizontal mo- 
tion. The horizontal detection screens (not shown) would 
then have to be vertically movable from close to the grat- 
ings to the far field, as described above. 



D. Dissociation-time entanglement 

One can also conceive schemes that go without grat- 
ings, by directly producing the modular momentum en- 
tangled states. Here we discuss a method based on 
dissociation-time entanglement [151 I16j . the sequential 
dissociation of a Feshbach molecule at two different times. 
A sequence of two dissociation pulses can generate a dis- 
sociation state where the atom pair is described by a co- 
herent superposition of an early and a late wave packet 
component associated with the two dissociation pulses. 
The two counter-propagating atoms are thus correlated 
in the dissociation times, see Fig.[3](b). 

The early and late wave packet components are spa- 
tially separated but propagate with equal velocities. 



They are thus described by a modular momentum en- 
tangled wavefunction with N = 2, similar to (13). The 



MME state is here realized in the longitudinal motion 
which separates the particles. The two dissociation-time 
components thus take the role of the slit components pre- 
pared by double slits, while the dispersive time evolution 
leads eventually to the overlap of the two wave packets. 

Position measurements in the overlap regions can be 
implemented by resonant photon scattering. The joint 
probability for the particle position then exhibits a non- 
local interference pattern. This completes the analogy 
with the Young double-slit experiment. To prove entan- 
glement one needs again a complementary measurement; 
in this case one must detect the positions of the atom 
pair at a time when the wave packets do not yet overlap 
to ensure that the particle pair is correlated in the early 
or late dissociation time. 



V. CONCLUSIONS 

We discussed a generic scheme to generalize the 
Young interference experiment for the case of two en- 
tangled particles, where nonlocal spatial interference is 
achieved by subjecting each particle to a grating struc- 
ture. The corresponding quantum state exhibits strongly 
non-Gaussian continuous variable entanglement, which 
can be revealed by a variance-based entanglement cri- 
terion. The latter is formulated in terms of modular 
variables, i.e. coordinates adapted to spatial interference 
phenomena. 

Experimentally, the entanglement detection is based 
on simple position measurements directly behind the 
gratings or in the far field; the modular variances are 
then calculated in a post-measurement analysis. We find 
that the entanglement detection scheme is quite robust 
against noisy EPR sources, coping even with substan- 
tial admixtures of classical correlations and incoherences. 
Moreover, while already double slits arrangements allow 
one to verify entanglement, one can improve the correla- 
tions by increasing the number of slits. 

We showed that a nonlocal Young test could be per- 
formed in a wide range of physical systems. Its demon- 
stration with material particles would be a striking 
achievement, demonstrating both the wave-particle du- 
ality and the non-locality of quantum mechanics at the 
same time. 
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